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Abstract 

The object of this paper is to investigate, classicahy and quantum mechanically, the relation 
existing between the position-dependent effective mass and damping-antidamping dynamics. The 
quantization of the equations of motion is carried out using the geometric interpretation of the 
motion, and we compare it with the one based on the ordering ambiguity scheme. Furthermore, 
we apply the obtained results to a Fermi gas of damped-antidamped particles, and we solve the 
Schrodinger equation for an exponentially increasing (decreasing) mass in the presence of the Morse 
potential. 
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I. INTRODUCTION 



In the course of the development of the quantum theory, Schrodinger's equation has 
played a central role It represents the basic tool for studying the behavior of many 
quantum systems ranging from nuclei j2| to sohd-state compounds js]. Very often, one 
deals with systems whose masses are constant in time and in space. Fortunately, this leads 
in many cases to exact results regarding the energy spectra and the eigenf unctions of the 
associated Hamiltonians, as is the case for, among others, the harmonic oscillator and the 
hydrogen atom. 

There exist, on the other hand, physical systems for which the mass exhibits the property 
of being time and/or position dependent. In particular, the motion of position-dependent 
effective mass (PDEM) has been the subject of a great deal of investigations, owing to its 



connection with the description of various physical systems such as quantum dots 



semiconductors |6l-ll3|. metal clusters 14|, ^He clusters 



ts 



15|, and quantum liquids IgI. De- 



pending on the intended application, many techniques have been used in order to study the 



PDEM motion |17l-l32j ; the standard approach consists in solving the Schrodinger equation, 
where one inevitably has to deal with the nontrivial ordering problem of the momentum 
operator and the mass in the kinetic energy part of the PDEM Hamiltonian 33|. 

In this paper, we explore the PDEM motion from a dynamical point of view, and we 
investigate its damping-antidamping character both classically and quantum mechanically. 
Generally speaking, the damping of classical particles is described by a frictional force that 



is proportional to the velocity (i.e., a force which is an odd function of the velocity) [34 1. 
In effect, the system under investigation is dissipative and is, quite often, non-Hamiltonian, 
making it difficult to quantize the equations of motion. This is the reason for which different 
approaches are adopted in the attempt to achieve a satisfactory quantum mechanical de- 



scription of damped systems 
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36| . Moreover, frictional forces that depend on the square 



of the velocity are known from hydrodynamics. As we shall see latter, the PDEM motion 
is tightly related to such forces, which may display antidamping behavior due to the fact 
that they are even functions of the velocity. The question that naturally arises is whether 
the corresponding equations of motion can be quantized using the standard technique. We 
shall be concerned with the answer of this question in the first part of this paper. 

The Fermi gas has found many applications in various fields of physics, like nuclear 
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physics js^], and solid-state physics jsl. Hence it is of interest to investigate the consequences 
of the spatial dependence of the mass on the properties of this system. In this spirit, we feel 
tempted to extend the analysis to the molecular domain, by studying the PDEM motion 
under the effect of the Morse potential ^8\. 

The paper is organized as follows. In Section Hi] we establish the relation between the 
PDEM and damping-antidamping motion. There, the quantization of the equations of mo- 
tion is carried out through the geometric interpretation of the dynamics, which enables us 
to find the stationary states corresponding to quasi-free particles for particular forms of the 
spatial dependence of the mass. Section [1111 is devoted to the application of the obtained 
results to a gas of damped-antidamped particles at zero and nonzero temperatures; the 
emphasis will be on the consequences on the pressure and the specific heat of the gas. In 
Section IIV| we solve the Schrodinger equation when the particle is endowed with exponen- 
tially decreasing (or increasing) mass, under the effect of the Morse potential, to determine 
the energy spectra and the eigenfunctions. The paper is ended with a summary. 

II. DAMPING-ANTIDAMPING INTERPRETATION OF THE PDEM MOTION 
A. Classical treatment 

In one dimension, the Hamiltonian corresponding to a position-dependent effective mass 
m{q) = moM(g), where M{q) is a dimensionless positive function of the generalized coordi- 
nate q, is given by 

Here V{q) is the potential, and p denotes the generalized canonical momentum conjugate to 
q, that is: 

{q, P} = 1, {q, q} = {P, P} = 0, (2) 

where the curly braces denote Poisson Brackets. Using the canonical equations of motion 

dH 

dp ' 



it can easily be shown that 



i^ = ^[V,m(?)]g'-V,V^(g), (4) 



where Vq = dq, and the dot designates the derivation with respect to time. Taking into 
account the fact that m(g) = qVqm{q), it immediately follows that 

Now multiplying both sides of the latter equation by the constant mo yields 



moq 



mo 
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VqM{q) q'-Vqi / -T7r^VqV{q)dq + Vo}, (6) 



2M{q) " '^^J^ " W M{q) 



q 



where Vq is a constant. Clearly, is the equation of motion corresponding to a particle of 
constant mass, mo, moving in one dimension, subjected to the effect of two forces; the first 
one being proportional to the square of the velocity: 

while the second one is a conservative force, derived from the effective potential 

Vesiq) = J + ^. (8) 

Integrating by parts, we can rewrite V^g as 

q 

Let us now take a look at equation ([7]). It is clear that is an even function of the 
generalized velocity q, meaning that the nature of the force is dictated by the explicit forms 
of M{q) and V{q). Suppose that the latter is such that the particle is constrained to move 
in some definite region of space. Consequently, is of frictional nature (the particle is 
decelerated) when g > and VgM(g) > or g < and VqM(g) < 0. On the contrary the 
particle gains acceleration (antidamping) if g > and VgM(g) < or g < and VgM(g) > 0. 
This may be better seen by noting that VqM(g) = M(g)/g. When V{q) is zero and M(g) 
monotonic, the particle (which we call quasi-free) is either damped or antidamped depending 
on its initial velocity. We thus come to the conclusion that the PDEM motion is equivalent 
to the damping-antidamping dynamics of a particle with constant mass subjected to a force 
proportional to the square of the velocity. 

In what follows we shall prove the converse. To this end, consider the following equation 
of motion 

moq = -(j){q)q^-VqU{q;a), (10) 



where U{q; a) is a potential depending on some parameter a. Our task here is to find the 
Hamiltonian from which equation ffTOj) is derived. First of all notice that a first integral 
for ([TU]) . which we denote by /, is a function of q and q, satisfying 

dl{q,q) _ -dl ^ ..dl _ Q 
dt dq dq 

This means that 



dl 1 



oq mo 



1 dl 

0(g)g' + V,[/(g;a)J— = 0. (12) 

Let C denote the family of curves in the (g, q) plane, for which dl /dq = holds. Hence along 
these curves we have that 

dq ~^ dq dq ^ ^ 
Comparing (fT2|) and (fT3|) . we get the following first-order ordinary differential equation 

^ = -—Mq - ^V,[/(g; a). (14) 
dg mo mog 

We seek a solution in the form 

g = /(g)e-^-^'^'^('^)'^^ (15) 

where /(g) is some function of g. By inserting the proposed solution into (fT4|) . we get after 
some algebra 

^!(! = _Ae4/.*(^)'^'?V„f/(g;«). (16) 
cig mo 

This can easily be integrated, and we get 

/(g)2 = _ A /" '^'•"^'"'VMiq; a)dq + C. (17) 
rriQ J 

Here C is a constant of integration characteristic of the curves C. Finally, the general solution 
of ([HD is 

g2 = Ce-^o-^o - ^ /" e^ *('^^''^V„t/(g; a)cig. (18) 

We solve for C to obtain 

C(g, g) = g^e^ + ^ L^J. "^("^'^ V,f/(g; a)dg. (19) 

<? 

The simplest choice we can make for the first integral of the equation of motion is 

/(g, q) = \mog{a)C{q, q) + K{a), (20) 
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where g{a) is a dimensionless positive function of the parameter a and where K{a) has the 
dimension of energy. 

According to j39-4^, the Lagrangian of the system is given in terms of the first integral 



by 



L = q I —^dq. 



A straightforward calculation yields 



L = -mog{a)q e™o 



VqU{q; a)dq — K{a.) 



(21) 



(22) 



It is now easy to show that the generalized momentum is equal to 

p = — — = mog[a)e"^o Jg^'-^' ^g. 
oq 

Consequently, the Hamiltonian corresponding to (|TOl) is given by 



H{q,p) = qp-L 

p2 



-— f Mq)dq , / \ I — f <t>{q)dq 



2mog{a 

Clearly, this is the Hamiltonian for a position-dependent mass 

m{q) = mo5'(a)e™o 
subjected to the effective potential 

U,f,{q;a) = g{a) j ^^'^^''^'''VgU{q;a)dq + K{a). 



(23) 



VgU{q;a)dq + K{a). (24) 



(25) 



(26) 



As an example consider the motion of the particle under the effect of the force —rjq"^ 
{4^{q) = f] > 0), in the presence of the Morse potential 



f/(g; a) = A(e~2ag _ 2e-"'?) 
with A,a > 0. From equation fl26|) . it follows that when a ^ rj/rriQ and a ^ 27] /rriQ, 

2'q 



(27) 



[/eff(g;a) = g{a)A 



exp{2(^-a)g} 2exp{(^-a)g} 



1-^ 

mo a 



2r? 



+K(a). 



(28) 



On the other hand if a = rj/mo, then 



C/eff(g; a) = 2g{a)A exp{2aq) — aq + K{a 



(29) 
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whereas for a = 2T]/mQ we find that 

Ues{q; a) = 2g{a)A exp{—2aq) + aq + K{a). 



(30) 



Obviously, g{a) and K{a) are not the same for all the cited circumstances; note, however, 
that while K{a) may arbitrarily be chosen (additive constant), g{a) should be reduced to 
unity when — t- 0. Even though, the choice one can make for g{a) still seems to be not 
unique. It should be stressed that, generally speaking, different Hamiltonians describing the 
same classical system may yield completely different quantum behaviors. 

To illustrate the importance of the inclusion of the function g we calculate the minimum 
of the effective potential Ues appearing in equation (128|) : 

1 2 



rrmm 



Ag{a) 



1- ^ 

moa 



moa _ 



-K{a). 



(31) 



We observe that the quantity Q 



blows up in the neighborhoods of 



1- — — 1-2- 
moa moa 

a = rj/rrio and a = 2ri/mo. Indeed, Q — )■ — oo when rj/{mQa) — t- 1+ or 2ri/{moa) — )■ 1~, 
whereas Q — t- +oo when ri/{moa) — ?■ 1^ or 2ri/{mQa) — t- 1+, which, obviously, is physically 
unacceptable. The only way to remove these singularities consists in choosing the following 
form for the function g{a): 



9{a) = g[a, 



1 _ ^)(i _ ^ 

mna ' ^ mna 



moa 



moa ' 



(32) 



where g[pi) is a regular positive dimensionless function of a. Therefore the effective potential 
reads 



t/eff(g;a) = g{a)A 



sgn( 1 



moa 



1 



2r? 

moa 



2sgn(l-^ 



moa 



V 
moa 



+ K(a). 



(33) 



As a result, the mass becomes a dependent. 



B. Quantization 



From a geometric point of view, the Hamiltonian ([T]) describes, strictly speaking, the 
motion of a particle with constant mass tuq in a one- dimensional space whose metric is 
defined by 

ds^ = M{q)dq\ (34) 



The corresponding Hamiltonian operator reads 



H = 



2mr 



(35) 



where A denotes the Laplacian operator. Prom the differential geometry we know that the 
Laplacian operator in the curvihnear coordinates with metric tensor g^^, is given by 



1 



E 



d 



det(g/,^)^. 



-1 d 



It follows that the time-dependent Schrodinger equation, in our case, takes the form 



(36) 



d 



d 



2mo ^JM{q) dq ^M{q) dq 



+ V{q) 



dt 



(37) 



In turn the wave function ■0(g, t) should be normalized according to 



\^{q,t)\'y/M{q)dq^l, 



(38) 



which actually arises from the fact that the line element is equal to ^/W{q)dq. On the other 
hand we can write: 



V dt ^ dt J 2y/^ 



\/m{q) 



m{q)^/'^ ■m{qy/'^ 



l V,m(g) _ _J^^ I V Ml) ^ 



2^/m{q) \y/m{q) 



ym(g) 



(39) 



where = dq, and = d^. Now by noting that the gradient in the q direction is 
grad = = M(g')~^/^Vg, the equation of local conservation of probability takes the form 



~dr 



with the probability current 



2im,o 



[r{q,t)Vqij{q,t) - ij{q,t)Vgr{q,t)]. 



(40) 



(41) 



Moreover, if we assume that m(g) does not depend explicitly on time, then we get another 
continuity equation in the usual Cartesian space: 

dp 



dt 



(42) 



where 



and where 



p{q,t) = mq,t)\'./M{qj = mq,t)M{qy/^\ 



(43) 



2imQ^M{q) 



h 



2imn 



r 



r 



M{qy/^ " \ M{qy/^ / M(g)V4 1 1 M(g)V4 



(44) 



Hence it is constructive to define a new wave function 



(45) 



M(g)V4- 

After some algebra, it can be shown that the latter wave function satisfies the Schrodinger 
equation 

- wi%) + ^^^^^ = '^W 



where the effective potential is 



V{q) = Viq) + 



]moM{q) 



h{V,M{q)f A,M{q) 
4M(g)2 M(g) 



(47) 



Case of a quasi-free particle in one dimension 

Suppose that the potential V{q) is identically equal to zero in the whole space, and that 
the position dependence of the mass of the particle is given by 



m{q) = mo e ""o 



(48) 



We shall look for the stationary states satisfying the time-independent Schrodinger equation 



Hi;iq) = Ei;iq), 

where > is the energy of the particle. Taking into account (1371) we get 

d'^ip 7] dip 2mo-E 



+ 



+ 



-(2r]/mo)q^ _ 



(49) 



(50) 



dq'^ tuq dq h? 

By the change of variable ^ = e"^''/"^"-''', the above equation reduces to the simple second- 
order differential equation 

d'^ilj 2mlE 
9 



which can easily be solved. Hence the general solution of (150|) is 



m = ^i(^)exp|y^e-(^/-)^|+C.(,)exp|-y^e-(^/-k^ 

= CMnq) + C2iv)i^^{q). (52) 

The meaning of the notations 'ip^iq) is clear from the above equation. We redefine the 
constants Ci and C2 in such a way that ip takes the form (assume r] > 0) 



m = C^i(^)exp^-M/^(l-e~^^)l+C'2(r/)expL/^(l-e-^«)^. (53) 



Now the limit of constant mass (r/ — t- 0) may be obtained by a Taylor expansion of the 
exponential function; this yields 

lim 7/;(g) = Ci e-''"^ + C2 (54) 



where k = ^Jlm^ElW' and where C\^2 = lini,,_>o (^1,2 (^)- It is clear that the above equation 
represents the usual expression of the wave function for a free particle. 

Consequently, we can interpret ip^{q) as the wave function of the particle when it is 
moving to the right, whereas ip~{q) describes the motion in the opposite direction. Indeed, 
it can easily be verified that the probability current density corresponding to ip"^ is equal to 



±y/2mQE/{h). Since Vgm(g) < 0, we conclude that ip' corresponds to the damping of the 
particle; on the contrary, describes the antidamping behavior of our quantum system. 

Now let us examine the transmission and the reflection of the particle from the rectangular 
potential wall U{q) = UQ6{q — a) where Uq > and 6{q) is the heaviside function. We further 
assume, for the sake of generality, that the mass of the particle is given by 

m{q) = moexp{-2(?7i6'(a - q) + r]29{q - a))q/mo}. (55) 

Note that the solution of the Schrodinger equation when the potential is U{q) = Uq = const 
can be obtained by simply replacing E in ( !52|) hj E — Uq. 

First we need to determine the conditions that should be satisfied by the wave function 
and its derivative at the discontinuity point a. To this end let us multiply both sides of 



equation (H9|l by ^y^I{q) and integrate between a — e and a + e, where e is a small number; 
we find 



lim 



2mo e-s>o 



dtp 



^M{q) dq 



dip 



^jM{q) dq 



ra+e 

lim / [E- U{q)] ^M{q)dq. (56) 

J a—e 
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The right-hand side of the latter equation is identically equal to zero since the quantity under 
the integral sign is bounded (we assume that M(g) is bounded in the interval [a — e, a + e]). 
Hence 



dip 



1 



dip 



(57) 



VM(g) dq J ^_ \,/M{g) dq 
which simply expresses the fact that the gradient of the wave function should be continuous 
at the discontinuity point. The continuity of the wave function itself follows from 

"a+e ^ 



= lim 

e-S>0 



..^^/M{^dq = iP{q)U+-iP{q)\a- 



(58) 



: y/M{q) dq 

Taking into account the above boundary conditions, together with the fact that the prob- 
ability current corresponding to ip^ is independent of the position q, we can show that the 
reflection coefficient has exactly the same form as in the constant mass case, namely: 



R 



1- Jl-t 



1 + 



1 _ Uo 
^ E 



(59) 



independent of rj. 

The Schrodinger equation we have used thus far was derived from the geometric inter- 
pretation of the PDEM motion. Nevertheless, it is worthwhile to point out that in dealing 
with the PDEM problems, one often writes the kinetic energy part of the Hamiltonian in 
the form 



T 



4mr 



-[M{qYPM{qYpM{qY + M{qYpM{qYpM{qY]. 



(60) 



Here /x, and k are von Roos' ambiguity parameters, which satisfy the condition ^ + v + n 



-1 



33| . Let us consider, for instance, the following Hamiltonian {V = 0) 



H 



1 



1 



-P- 



(61) 



2mo M(g)' 

obtained from (1601) by setting u = k, = 0. It can easily be verified that the Schrodinger 
equation with this Hamiltonian yields 



dq'^ mo dq 



V^ = o, 



(62) 



which differs from equation fl50|) by the factor 2. Now by making the change of variable 



(63) 
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To eliminate the minus sign in front of the first derivative we further assume that = C^P- 
We then find that the function (p satisfies 

e^" + e0' + ie - 1)0 = 0. (64) 



This is the well-known differential equation for the Bessel functions of order 1 43(| ; thus the 
solution of f l63|) is: 

m = C,{v)Ji{0 + C2{v)Yi{0, (65) 
where Ji and Yi denote, respectively, the Bessel functions of the first and second kinds. It 
follows that the general solution of equation (l62l) is given by 



"(66) 

Besides the fact that this wave function describes a completely different quantum behavior 
as compared to the one obtained above in (|52l) . here the smooth transition to the constant 
mass case cannot be made. Notice that this transition was possible in flS21) because of 
the multiplicative property of the exponential function, which does not hold for the Bessel 
functions. In particular, the explicit form of the probability current density corresponding 
to fl66|) turns out to be quite involved. 



III. FERMI GAS OF DAMPED-ANTIDAMPED PARTICLES AND EFFECTIVE 
MASS TENSOR 

In this section we generalize the concept of the Fermi gas to the case where the particles, 
which we suppose identical with mass m (we omit for now the subscript 0), are individually 
subjected to the effect of three forces; each of these is proportional to the square of the 
component of the mechanical velocity along the corresponding coordinates axis of the three- 
dimensional configuration space, that is, 

Fq = VVq, <l = x,y,z. (67) 

Using the results of the previous section, we find that such a system admits the Hamil- 
tonian 

N / 2 2 2 \ 

^l2m 2m 2m J' ^ ^ 
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where is the number of particles. Accordingly, we may regard the g system of 

particles in a medium where each particle is characterized by the diagonal effective mass 
tensor diag(me~^'''^/™, me"^™/*", me"^''^/™). 

We assume that the gas is in thermal equilibrium at temperature ksT = 1/(3 {ks denotes 
Boltzmann's constant) and is completely contained in the cube —L < x,y, z < L. This is 
the only acceptable form of the intervals, otherwise the symmetry with respect to the change 
of the sign of the force will be violated. Under the latter assumptions the classical partition 
function reads 

3N ^ oo 



k=l 



/27rm\3JV/2 




(69) 



Clearly, V{ri) represents the volume in the 3D curvilinear space with the metric tensor 



diag(e 



-rix/m, p-r)y/m 



e '?^/™)^ which justifies once more the use of the geometric interpretation 



of the motion. The integral in (l69l) can easily be carried out, yielding (for t] > 

3N 



z 



1 /27rm^37V/2 



sinh {r]L / m 



(70) 



Hence, by Stirling's formula, the free energy of the system can be written as 



—kBTlnZ 
-kBTNi ?lnT + ln 



{2m /r])^ sinh'^ (r]L/m 



+ In 



{2'KmkBf^^t 



(71) 



Knowing the explicit form of the free energy, one can calculate the diverse thermodynamical 
quantities characterizing the gas. In particular, the pressure turns out to be equal to 



with 



46| 



P 



V{v) 



/dF\ _ NkeT 
XdVJTN ~ V{ri) ' 



sinh i^L/m 



Thus 



PV{r]) = NkBT, 



(72) 



(73) 



(74) 
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implying that at constant temperature, the pressure decreases as t] increases (see discussion 
bellow). The chemical potential, on the other hand, reads 



2 \/ 2-11 k bTti^T sinh yrjL / m 
In a similar way we find that the internal energy of the gas is 



(75) 



which does not depend on i]. As a direct consequence, we conclude that the specific heat 
is equal to Cy = ^Nks- Hence, although the pressure is inversely proportional to 1], the 
classical analysis reveals no dependence of the internal energy and the specific heat on the 
applied force, in contrast to what one might expect. As we shall see bellow the quantum 
mechanical treatment of the problem does resolve this issue. 

To achieve this aim, we first look for the discrete energy levels of the gas, assuming that 
the potential is zero within the region of space containing the particles and is infinite outside 
it. The wave function of every particle of the gas is simply given by the product of the wave 
functions associated with the one-dimensional motion in each coordinates axis, that is 

^(a:, z) = ijix)ip{y)^/j{z), (77) 

where 



m = C^iv)ewi^^^^e~^^/-^''yCMew\^^^^ 1 (78) 

Since this should vanish at the edges of the interval —L < q < L, it follows that 



which yields the quantized energy levels (g = x, y, z) 

En, = ( Si 5"' , , , r., = 1,2,3,... (80) 



Hence the eigenenergies corresponding to fl77|) read 



En^uyu. = b ^« + ^D- (81) 

8m^sinh irjL/m 
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The wave function (ITHj) . in turn, takes the form (we drop the subscript q) 



^n(g) = C(")(r7)<|exp| -^J^^^e-^"/'")'' 1-expi 



U-Vq/m_2Q-vL/r, 



^2^2 ) \ ]l ^2^2 

The constant C^"'\r]) can be determined from the normahzation condition 



|V^n(g)re-^^/"^cig = l. 



After some algebra, one obtains 



~r)L/m 



cos[2k{En)e 



sin I 2k{En)e 



-•qL/m 



-■qL/m 



cos(2k{En)e 



TjL/m 



2k{E„ 

-1/2 



cosi2k{En)e 



sin I 2k{En)e 



r)L/m 



— I 2 sinh (nL I 'm\ — sin 

2m I V ' } 2k{En) 



-1/2 



4:k{En) smh(r]L/m 



with k{En) = ^2m?En/h^ri'^. 

For completeness we also display the quantity 

L 

Vn{L) = [ \MQ)\'dq 



'^\C^^)(n)\^!^2L + sin(2A;(E„)e-''^/"^) 



X 



X 



Si (^2fc(E„)e-'^^/™) - Si (^2fc(E„)e''^/'")) + cos (2fc(^„)e 



Ci(2A;(E„)e-^^/™) - Ci(^2fc(E„)e^^/" 



-rjL/m 



where 



Si(a;) 



Ci(x) 



sin(0 



cos(0 
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are, respectively, the sine and the cosine integral functions [43|. In figured] we display the 
quantity p„(g) = \tpn\'^ ^/ M (q) as a function of q for some values of rj, L and n. We see that 
the probability of finding the particle is maximum in the region of space where the mass 
is large (see the discussion bellow). The same observation holds for |V^„p, as indicated in 
figure El 
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■3- 2.0 
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0.0 




ii=4.L=].);= 




FIG. 1: (Color online) The probability density Pniq) = \'^n\'^ \/ M {q) as a function of position for 
some values of L, ry and n with m = h = 1. 



At zero temperature all the levels are occupied up to the Fermi level, ep'-, hence if we 
assume that the former are quasi-continuous |3|, then the density of states is equal to 



dN 



sinh^ (^r]L /rnj VEdE, 



where g denotes the degeneracy. Thus the total number of particles reads 



sinh^ (^r]L/m^ y/EdE, 



from which we find that 



2^2 



hi] 



1 



67r2A^ 
9 



sinh^ ( rfL / m 



(88) 



(89) 



(90) 
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FIG. 2: (Color online) I'i/'nl^)!^ as a function of position for some values of L, rj and n with 
m = h = 1. 

This equaUty reveals that the Fermi energy gets smaller as r] increases, in accordance with 
relation fIST]) . Now the total energy of the gas can be calculated using 



which yields the familiar expression 



Etot = j EdN, 





Etot = -Nep. 
5 



(91) 



(92) 



The fact that the total energy is still proportional to the Fermi energy explains the reason 
for which the pressure of the gas decreases due to the increase of rj, as was found earlier via 
the classical partition function. 

Actually, using the results of the previous section, together with pure classical arguments, 
we can explain this behavior as follows: For every particle the conservation law {q = x, y, z) 



(93) 



holds at any moment in time. Hence, for 77 > and q close to L, the velocity of the particle 
is quite large which means that it quickly leaves this region just after hitting the wall at 
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q = L. On the contrary, the velocity is quite small in the neighborhood of —L: the particle is 
thus slowed down in this region, taking larger times to reach the wall at g = —L. Although 
the carried momentum is large near q = L, after sufficiently long time, there will be much 
more particles close to —L, as clearly indicated in figure [T| hence, as a whole, the pressure 
of the gas drops down. 

At T 7^ 0, the number of particles is given by 



N 



g^9/2^ 



sinh'^ ( rjL / m 



EdE 



l + exp[{E - fi)/kBTy 



(94) 



where yU denotes the chemical potential. By introducing the notation z = explfi/ksT), and 
making the change of variable ^ = exp{—E/kBT), we can write 



N 



sinh irjL/m 



Integrating by parts yields 

N = _^-^'''9{ksTf'' sinh3f,L/m)Li3 (-e^/^^^. 

^^3/2^3^3 V W ) 

Here, Li„(a;) is the polylogarithm function, defined by 



(95) 



(96) 



\n — 1 



Li„(x) 



Tin - 2) 



ln"-'(2/)ln(l-a;y)— . 

y 



(97) 



Then using equation fl90|) we find that the chemical potential is given in terms of the Fermi 
energy by 



/i = kBT ln< — Lig ^ 



3/2- 



(98) 



where the superscript designates the inverse function. Similarly, the total energy of the gas 
at nonzero temperature is found to be equal to 

6m9/2^(fcBT)5/2 



E,. 



V^7r3/2;i3^3 



sinh^ {r]L/m^ Li s (-e'^/^^^) . 



9^NkBT ( keT 



3/2 



Lis(-e'^/^^^). 



Using the property j44| 



:^Li„(x) = 

ax X 



(99) 



(100) 



18 



we find that the specific heat of the gas is 



V 



dT Jv 



3/2 



v 



5k 



Li,(-e^ABn + i^Li3(-e^/'=-n 



For T ~ we may use the the approximations 45 1 



1-^ 



12 \ ep 



tot 



\Nep 
5 



1 + 



StT^ / knT 



12 V ep 



This yields 



Cv ~ ^'K^NkB-, 



which implies that the specific heat of the gas increases with r]. 



IV. MOTION IN THE MORSE POTENTIAL 



(101) 



(102) 
(103) 

(104) 



In this section we propose to solve the time-independent Schrodinger equation 
ri^ 1 d 1 d 



2mo ^yM{q) dq ^M{q) dq 



U{q;a) - E 



for the Morse potential 



(105) 



U{q;a) = A{e 



-2aq 



2e-"^), 



A, a>0. 



We shall consider two special cases, namely, M(g) = exp{±2riq/mo}, with 

rj/rrio = a. 



(106) 



(107) 



This choice leads to exact analytical solutions as we shall see bellow. Physically, the latter 
condition implies that the effective distance constant of the force ([7]) is equal to the effective 
range (1/a in lenght unit) of the Morse potential. 

Consider first the case M(g) = exp{2?7g/mo}, and let us look for the discrete spectrum 
corresponding to negative values of E. By virtue of equation f ll05p . we can write 



d'^ip 1] dip 
dq^ niQ dq 



+ 



2moA 



= 0. 



(108) 
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Making the change of variable 



(109) 



equation fllOSp becomes 

ip + 



1 
4 



2mlA \ 1 



2mlA \ 1 



^ = 0. 



(110) 



We require that the function vanishes as g — > ±oo; in terms of the variable ^, we have to 
study the behavior of ip at infinity and in the neighborhood of zero. One can show that 



if) r-u Q 



(l+Vl+4s2)/2 



-«/2 



as ^ — )■ 0, 
as ^ — 7- oo, 



(111) 
(112) 



where we have introduced the quantity 



(113) 



Hence it is natural to seek a solution in the form 



(114) 



where is a new function of ^. By direct substitution into flllOp . it can be verified that 



^0" + (1 + Vl + 4s2 - 00' + [r - (1 + Vl + 4s2)/2]0 = 0, 



with 



2mlA 



(115) 



(116) 



Equation (IllSp is nothing but the differential equation associated with the confluent hyper- 
geometric function 43|, |44 1 

= F(-n,2« + l, 0, (117) 



where 



K = ^Vl + 4s2, 



n = T — K. 



(118) 



It follows that the unnormalized eigenfunctions take the form 



^n(0 = e-n^^+^)/^F(-n,2«: + l, 0- 



(119) 
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The required conditions on the behavior of the wave function at ±00 imphes that only 
positive integer values of n are allowed, whence 



— Er, 



(n+|) + iJl + ^^ 



2 ' 



n = n. 



mm; ^mm 



+ 1;... 



(120) 



2^ ' 2 V " ' h^v'^ 

Here, nmin denotes the minimum positive integer number for which —E < A. Thus, although 
the spectrum of negative eigenvalues corresponding to the Morse potential consists of a finite 
number of levels for constant masses, the chosen form of the position dependence of the mass 
we have adopted above renders the number of energy levels infinite. 
Now if we assume that M{q) = exp{— 2?7g/mo}, then we have that 



dPifj T] dtp 
dq^ rriQ dq 



^ = 0. 



;i2i) 



By the change of variable 



we get 



(122) 



(123) 



2171^/1]'^ 

which is the differential equation for a parabolic cylinder function. Let us further introduce 
the notations 

o 2Ar]'^ , , mn 



rriQ T] 
Hence with = 1 — ^, equation fll23p takes the form 

d^ 

— \- 

2 



M 



-,2 ' 



£ = E + A. 



(124) 



(125) 



2Md^^ 

Clearly, the above procedure enables us to map the whole dynamics into a simple harmonic 
oscillator problem in the interval ] — 00, 1]. It follows that the solution of fll2ip is 



with the eigenvalues 



ml A 
2hhf 



E„ 



1 



-riq/mo 



prr4A 



h^|'^^{n + i) 



A. 



■rri'A 



(126) 



(127) 



In the above, C„ is a constant and Hn{x) denotes Hermite's polynomial of degree n; the 
allowed values of n are those which satisfy the conditions 



-Er,. < A. 



;i28) 
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The former equation can be solved either geometrically or numerically. Nevertheless, using 
the property 

i^2m+i(0) = 0, m = 0,l,2,---, (129) 

we conclude that for sufficiently large values of r], only odd values of n should be taken into 
account, provided —En < A. 

V. SUMMARY 

Summing up, we have shown that the motion of position-dependent effective mass can 
be regarded as a damping-dantidamping dynamics of a constant mass under the effect of 
an effective potential, and vice versa. Care has to be taken when calculating the effective 
potentials since certain singularities may appear, which affects, in turn, the mass of the 
particle. We have quantized the equations of motions starting from a geometric interpre- 
tation of the motion; we found that it allows for a smooth transition to the constant-mass 
case, which does not hold for the momentum and mass ordering method. We have applied 
the obtained results to a Fermi gas of damped- ant idamped particles. The classical and the 
quantum treatment of this system reveals that the pressure and the specific heat depend 
on the applied forces. Finally, by solving the Schrodinger equation in the presence of the 
Morse potential, we deduced the energy spectra and the eigenfunctions of the bound states 
in the case of a particle with exponentially decreasing or increasing mass. We found that 
the number of energy levels may become infinite for particular spatial dependence of the 
mass. 
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